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Abstract 


This study addresses the inverse issue of identifying the space-dependent 
heat source of the heat equation, which is stated using the optimal con- 
trol framework. For the numerical solution of this class of problems, an 
approach based on shifted Legendre polynomials and the associated oper- 
ational matrix is presented. The approach turns the primary problem into 
the solution of a system of nonlinear algebraic equations. To do this, the 
temperature and heat source variables are enlarged in terms of the shifted 
Legendre polynomials with unknown coefficients employed in the objective 
function, inverse problem, and initial and Neumann boundary conditions. 
When paired with their operational matrix, these basis functions provide 
a quadratic optimization problem with linear constraints, which is then 
solved using the Lagrange multipliers approach. To assess the method’s 
efficacy and precision, two examples are provided. 
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1 Introduction 


Many inverse problems for the heat equation are applied to many fields of 
physics and engineering, such as acoustics [19], medical imaging [10], signal 
processing [35], optic [6], and radar [7]. There are approximately five major 
classes of inverse heat diffusion equation problems. 

i) The problem of reverse time or conducting heat backward from the known 
last-minute distribution determines the initial temperature distribution. 

ii) Inverse heat conduction is the detection of temperature or temperature 
flux at one inaccessible boundary beyond the data available in the other case 
that is accessible. 

iti) Identify coefficients of over-posed data at the boundaries. 

iv) Determining the shape of unknown boundaries or cracks inside the heat 
conduction body. 

v) The identification of the heat source [18, 5]. 

The heat equation, in this research, treats the heat source as an uncertainty. 
Applications in the real world where these difficulties are useful include cre- 
ating the end state of melting, and freezing processes and determining the 
contaminating source’s intensity. Methods such as the generalized finite dif- 
ference scheme [12], the radial basis function method [28], the sparse reg- 
ularization approach [25], the meshless generalized finite difference scheme 
[13], the mollification regularization scheme [36], and the reproducing ker- 
nel space scheme [33], have all been applied to the solution of inverse heat 
source problems. In this paper, we propose a novel numerical method for ob- 
taining the source parameter (or control parameter) in parabolic equations. 
Iterative methods and a variational approach have recently been proposed 
to numerically solve this problem [21]. These methods are computationally 
expensive because they solve a direct problem at each iteration. Tikhonov 
regularization is proposed in [37] as a stable optimal control solution to the 
inverse heat source problem. Parameter identification for a nonlinear heat 
equation in the 2D and 3D space-time domains was solved by Lin and Liu us- 
ing homogenization functions as the basis [24]. The authors of [29] proposed 
a perfect method to investigate inverse heat source problems in function- 
ally graded materials using the homogenization function. Due to the given 
conditions, a homogenization function for the boundary value problem is con- 
ceived, and a family of homogenization functions is further derived. Djennadi 
et al. [11] employed the expansion method and the overdetermination con- 
dition to solve the inverse source fractional diffusion problem that contains 
the Atangana—Baleanu—Caputo fractional derivative. In [20], for the stable 
reconstruction of the heat source in the parabolic heat equation, an itera- 
tive variable conjugate gradient algorithm is proposed based on a sequence 
of direct problems that are solved using the boundary element method of 
each iteration step. The gradient descent along with the finite difference 
method to find the solution nonlinear inverse heat transfer problem in [4]. 
Ciofalo [9] proposed using finite volume discretization to get a solution for 
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an inverse heat conduction problem. In which, with the assumption that 
the thermal boundary conditions in other walls are known, the steady state 
distribution of the displacement heat transfer coefficient on one slab wall is 
reconstructed from the temperature distribution in the plate embedded in 
the slab. With the increasing use of machine learning techniques, includ- 
ing neural networks, the use of these techniques in solving inverse problems 
has also attracted the attention of many researchers. Li and Hu [23] used 
a multi-layer neural network to solve the Cauchy inverse problem. Phisyc- 
informed neural network models are one of the powerful methods in deep 
learning. Authors in [27] applied it to solve a class of inverse problems re- 
lated to partial differential equations (PDEs). The authors of [15] proposed 
a new method for solving large-scale inverse problems based on Bayesian in- 
ference, Markov chain Monte Carlo approach, and derivative-free algorithms. 
Bondarenko [8] presented a finite-difference-based method to investigate the 
discrete systems of the inverse of the Sturm—Liouville problem. Huntul [16] 
used the Tikhonov regularization and the nonlinear optimization for the first 
time in the third-order pseudo-parabolic equation with initial and nonlocal 
periodic boundary conditions derived from nonlocal integral observation for 
the inverse space-dependent heat problem. Huntul [17] recovered a source in 
a high-order pseudo parabolic equation using cubic spline functions. In [14], 
authors solved the two-dimensional inverse time-fractional diffusion prob- 
lem with nonlocal boundary conditions using a-polynomials, collocation, and 
least squares methods. They calculated time using the Z; method. Wen, 
Liu, and Wang [34] used the Fourier approach to find the source term and 
starting data in the time-fractional diffusion equation . Abbaszadeh and 
Dehghan [1] considered the inverse tempered fractional diffusion equation. 
They used Crank—Nicolson temporal discretization, a modified element-free 
Galerkin method, and a meshless method to solve the inverse problem. 


In this research, we provide a numerical solution for solving the inverse 
heat source issue in an optimal control setting by using orthogonal polyno- 
mials. This piece is an attempt to provide a fresh strategy for addressing 
the issue of the mysterious heat source. The optimal control issue is re- 
duced to a set of algebraic equations in the suggested approach [26, 32, 30]. 
This is achieved by approximating the temperature y and the heat source 
f in Py (see (1)) with the help of shifted Legendre polynomials (SLPs) and 
their operational matrix with unknown coefficients. By substituting these 
approximations for the objective function in the inverse problem, we are 
able to determine not only the unknown coefficients but also the initial and 
boundary conditions. To conclude, we utilize Lagrange multipliers to connect 
the algebraic equation produced from the objective function to the algebraic 
equations derived from the inverse system and the starting and boundary 
conditions. Then, we can use the constrained extremum method to solve the 
resulting algebraic system of equations to find the best solution. The authors 
of [3] investigated the inverse heat equation problem with variable boundary 
conditions using a weak solution strategy. The Legendre spectral collocation 
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method was used to solve a fractional inverse heat conduction problem in 
[2], where both the temperature function and the boundary heat fluxes were 
unknown. Following the introduction, the article will focus on five primary 
sections that together address this inverse problem. In Section 2, we present 
the optimal control issue and the inverse heat source problem. In Section 3, 
we describe the SLPs and their characteristics. The problem is resolved in 
Section 4. In Section 5, we provide numerical examples that demonstrate the 
effectiveness and precision of the suggested approach. The last part explains 
the results. 


2 Problem statement 


Suppose the following inverse problem: 
Let us suppose 9 := (0,1)x(0,T), T > 1, one is going to find the temperature 
z and the heat source f that satisfy (1); that is, 


24(x,t) — Za (x,t) = f(x), (x,t) E€ 0, 
Pi: ¢ 2(z,0) = v(z), ee (0,1), (1) 


2_(0, t) = go(t), 2_(1,¢) = gilt), te (0,2). 


The second-order parabolic equation Py, with sufficiently smooth functions 
v(x) (the initial condition), (Neumann boundary conditions) go(t) and gi(t), 
forms the governing equations. 


Assume that the desired function measured data h,(a) (desired function) 
and the actual data z(x,T) := h(x) meet the following relation: 


I[h(a) — he() [22501] S & (2) 


where € is the known noise level and the norm || - || ,2/0,1) of a function z(2) 
is determined by 


Ile = (f° Peae)’, 


In the following part, we convert the problem Py; into an optimal control 
problem of Py and solve it using the suggested approach. The following is a 
consideration of the optimal control problem: 


: 1 o 
Pu: ymin J(2,f) = 5lle(e,T) — helo) aoa +SIVE laps @) 
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where F,¢ has defined the set of admissible controls of the objective function 
J as 


Fua = {f(x):0< as f <b,Vf € L[0, 1}, (4) 


with the constant bounds, a,b € R. Moreover, z(x,t) is the solution of (1) 

for a given heat source f(x) € Fyua, and o is the regularization parameter. 

For noisy data h,(a), the purpose of the optimal control problem is to find 

functions f(a) and z(a,t) that minimize the objective function Py and satisfy 
I- 


3 Shifted Legendre Polynomials (SLPs) 


The orthogonal polynomials with regard to the weight function W(x) = 1 
on [—1, 1] are known as Legendre polynomials of degree m and are denoted 
by Lm(x) (m = 0,1,...). The following recurrence formula can be used to 
create these polynomials: 


2m+1 m 
Lim = De, — —— 
(x) mea (x) aR 


Lini({2), m=1,2,..., (5) 


where Lo(a%) = 1 and Li(x%) = x. The well-known SLPs in [0,1] can be 
created by changing the variable x = 2t—1, which is expressed as £L,,(t) (m = 
0,1,2,...,), by 


(2m + 1)(2¢ — 1) m 
m+1 Em(t) 


Pa Lma(t), m=1,2,..., (6) 


where £o(t) = 1 and £,(t) = 2t — 1. The explicit formula of the SLPs is as 
follows [31]: 


Lim (t) = s binit’, (7) 
i=0 
where £,,(0) = (—1)™, £Lm(1) = 1, and 
. (m+i)! 
bee = (—1™ oe 8 
a) (m — 2)! (2!) .) 


The orthogonality condition of the SLPs with respect to the weight function 
w(t) = 1 is given by 


/ * Lpa(t)Ln(t)dt = indian, (9) 
0 
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where émn is Kronecker’s delta function and hy, = . Any given func- 


1 
2m+1 
tion z(t) € L?[0, 1] can be represented in (n + 1) terms of the SLPs as 


n 


2(t) ~ S zLi(t) 2 ZT, (t), (10) 
1=0 
where 
Z = [20 21 Zn)" 
&,,(t) £ [Lo(t) Li(t) ... Lalt)], (11) 
and 
= i 2(t)L;(t)dt, i=0,1,...,n. 
i JO 


In a similar way, a two-variable function z(x,t) € L?(@) can be expanded by 
the SLPs as 


2(a,t) ~ S\S> zijLi(a)L£j(t) 2 &2,(x)Z®,,(t), (12) 


i=0 j7=0 


where Z = [z;,] is the matrix of coefficients with dimensions (m+ 1) x (n+1) 
whose entries are unknown and obtained from the following equation: 


1 1 1 
ee | : 2(m,#)L;(a)L£, (deat, (13) 


for i = 0,1,...,m, 7 = 0,1,...,n. Suppose that ©®,,(¢) is the vector intro- 
duced in (11). Then the derivative of this vector is as follows: [31] 


d®,,(t) (1) 
—— =D; ’®,(t), 14 
where pW) = [ds] is called the derivative operational matrix of SLPs of 


(n + 1)-order, whose structure is as follows: 


k =1,3,...,n if n odd, 


2(2j+1), g=i-k 
(24j7+1), j=t ieee if n even, 


(1) _ 
a) = (15) 


0 otherwise. 
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Remark 1. Generally, the r-derivative operational matrix of SLPs of ®,,(t) 
can be given by [31] 
d’®,, (€) 
dt” 


= DIS, (t), (16) 


in which D\”) is obtained by r times multiplying D{” in itself. 


4 Convergence analysis 


In this section, the convergence analysis of SLPs expansion in two dimensions 
is investigated. 


Theorem 1. Suppose that z : 0 —> R is (n +m +1) times continuously 
differentiable. If ®7 (x)Z®,,(t) is a unique best approximation of z, then the 
following inequality holds: 


AVT(n+m +2) 
r(nt+m+1l—ryP 


|| (a, t) — F(x) Z®n(t)||L2(@) < (17) 
where 


ontmt+rl 


A= max{| x,t)| i=0,1,....n+m+1}, 


Oxr+m+1—t6 te 2( 
r= max{Ten Per", i=0,1,...,2(¢n+m-+1)}. 


Proof. Maclaurin’s expansion for z(«,t) reads as follows: 


1 O Op Son 
2(a, t) = p(a, t) (n+ m+ Dae + ta) a +1 2(€o2, £0t), £0 € (0, 1), 
where 

n+m k 
1, O O 
p(x,t) = » ries +ta) z(0, 0) 

Thus 

1 6) O 


Got m tI ae thay) kor Sot)], 0 € (0,1). 


|2(x, t) — p(x, t)| = | 


On the other hand, since ®7 (x)Z®,(t) is the best approximation of z(z, t), 
we obtain 


Il z— 82.29, |)? ,2¢6) Sill 2— PI? 226): 
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By the definition of the L?-norm and expand («2 +t#)"*™*!, we have 
l| 2(a, t) — p(x, t) IIZ2(e) 


T pl 
=[ ar 1 Ge 42. yntmt1 (eon, €ot)|2dardt 


m+1)!> da ot 
T n+m+1 
-/ fi 1 S- n+tm-+1 grtmt+1—izi 
@ so (emt a a 
gontm+1 


—.—.— — 2(€02, fot 2 dadt 
Bg ttt opi ( 0 £0 ) 
yntm+1 


A2 T _ 
Satie ff USE (Cnt aroma 
(n+m+1)? Jo Jo i 


where 


ea) =max{ ("1 *1)s¢=0,1,....n4¢m41}. 


To find the upper bound for the above expression, we calculate the following 
terms: 


T pl ne Plt2i 
nrm —14 t ‘d. dt = 
7 fe PTT + Di) (Qn + 2m — 2 + 3)’ 


7=0,1,...,n+m-+1, 
a i (Qn42m41—1) 4i41 ae 
nm 1 a dxdt = 
[ fe ve (24 (Qn + 2m —i +2)’ 
7=0,1,...,n+m, 


T 1 (244) ease : T2nt2m—i+1 
42n42m-i dade = = 071 
| | 7 "B+ )QntIm—it]ly |” 


T 1 ; : T2nt2m—i+2 
| / gilti)pant2m—itlogde = i=0. 
0 0 (2 7 


Therefore 


A2 T al ce 
| z — Pll’ 12,6) < 1p p | | [ gener dat 
r?(nt+m+1—r)? Jo Jo =i 
2 


n+m+1 Tit+2i 


: 3 
r?(n+m+1—r)!? = (1 + 22)(2n + 2m — 21 + 3) 
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n+tm Tesi 
» (2+ 4)(2n + 2m —i+2) _— 
1 ant+2m—i+1 Pant+2m—i+2 
* » + (Qn+2m+1) XIn+Im+ | 
Or 


< 
~ r2#(nt+-m+1—r)!? 


x [(n+mt2)+(n+m41)4+---+24]] 
Aer 
< 
~ r2(nt+-m+1-—r) 


g(a t+m + 2)’, 


which is the desired result. 


5 Description of the presented method 


Now in this section we will use numerical methods to address the problem 
raised in (1) and (3). We will use numerical methods to address the problem 
raised in (1) and (3) in this section. SLPs approximate the temperature and 
heat source for this purpose as follows: 


2(z,t) ~ On (x)? ZO, (t), (18) 
FT@,(z), (19) 


8 
~~ 
I2 


where Z and F are the following unknown matrices of coefficients with di- 
mensions (m+ 1) x (n+1) and (m+1) x 1, respectively, while ®,,,(x) and 
®,,(t) in (11) have been expressed: 


Z00 701 --- 20n f 
210 #11 +++ Zin h 
Z= ; = . (20) 
zm0 ant 2mn Jm 
Set 
P(a,t) = (Lo(x)Lo(t), --., Lm(z)Lo(t) | ++» | Lo(w)£n(t), --., Lm(x)Ln(t)]. (21) 
According to (21), we can express (18) as 
z(a,t) ~ ®,,(x)? ZG,,(t) = P(x, t) vec(Z), (22) 


where 


T 
vec(Z) = [zo00,---;2mo |---| Zons--+) Zmn 
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From (14), (22), and Remark 1, the result will be as follows: 


2,(x,t) ~ Bp, (2)? DY? ZO, (t) = P(x, t)\(Ing1 ® DY) vec(Z), (23) 


Znq (t,t) & ®m (x) DPT Z&n(t) = P(a,t)(Ing1 @ DP”) vee(Z), (24) 


x 


x(a, t) © Om (x)? ZD\S, (t) = P(a,t)(D.?” @ Imai) vec(Z), (25) 
so that Im+41 and I,41 are identity matrices of orders m+ 1 and n+ 1, 
respectively. Additionally, © refers to the Kronecker product [22]. Now, 


(19), (24), and (25) are substituted into the first subequation of (1), and the 
result is 


K(a, t) vec(Z) — F7®,,(x) = 0, (26) 
in which 
K(a,t) £ P(e, t) [(DP* @ Im) — Unt @ DO)]. 


Thus, as to (22) and (23) and with regards to initial and Neumann boundary 
conditions in (1), we have 


P(x, 0) vec(Z) = v(a), 
P(0,t)(In41 ® Dz") vee(Z) = go(t), (27) 
P(1,t)UInt1 ® DS?) vec(Z) = gi(t). 
We follow the suggested procedure by constructing an (m-+1) x (n+1) alge- 
braic system of equations. For this reason, we derive the following equations 


from (26) and (27): 


K(&, nj) vec(Z) — F? ®m(&) = 0, 1=2,...,m, J =2,...,n+1, 
P(&,0) vec(Z) = v(&) i=1,...,.m+1, 


P(0,nj)(In41 ® DE) vee(Z) = go(nj), § =2,.-..n +1, 


P(1,1j)Ing1 @ DY?) vec(Z) = gi(nj), §=2,---.2 #1, 
(28) 
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where a collocation scheme is defined by evaluating the outcome at the points 
(€;,7;) in (28). We employ the shifted Legendre-Gauss—Lobatto nodes €; (1 < 
i <m-+1) and the shifted Legendre roots n; (1 < 7 <n+1) of £,(t) to find 
suitable collocation points. It is possible to write (28) as follows: 


M vec(Z) —-NF = Q, (29) 
in which 
T(2:n4+1,:) @ &(2:m,:) S@X(2:m,:) 
_ T(1,:)@¥(1:m+1,:) 7 0 
ae T(2:n+1,:) ® X(1,:) : Me 0 : 
T(2:n4+1,:)@X(m+1,:) 0 
T 
Fe Fopiis5 Fel Opicsg 0 h[OpicegO)) 5 
—_—— —_——— 
m+1 m+1 
re Ooo 
Q=[0,...,0]-+-10,..-,01o(G1),---s¥(Ema1) 
ee SS 
m-1 m-1 
7 
|90(72),--++90(tm+1)|91(M2),---.91(Mm-+41)]| 
where F and Q are (m+ 1)(n + 1)-order vectors and 
10... 0 
ga 00... 0 
00 0 0 nx(n+1) 
Lo(m) £1(m) —--»  L£n(m) 
Lo(n2) £1(n2) L£n(N2) 
T= ; 
Lo(mn+1) Ly (itt) oo La (Ynat) (n4+1)x(n41) 
Lo(&1) £1(&1) L£m(E1) 
. Lo(&2) £1 (€2) £m (€2) 
X= ; 
Lo(Em+1) £i (Gms) see Plead) (m+1)x(m+1) 


Next, we approximate Py by SLPs. First, we approximate the desired func- 
tion h_-(x) as 
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h-(a) ~ H?®,,(a). (30) 
From (14) and (19), we have 
Vi (2) ~ FT DM 6,, (2). (31) 
Inserting (18), (30), and (31) into (3) yields that 
J(z,f) & Imn(Z, F) 
= sf @m(oy" 2000) — Om (x)7 H)(®m (x)? ZOn(T) — Om(x)? HH)? dx 


oO 


al 
a al (FT DY}, (2))(FT DY Bm (2)? de. 
0 


The value So (m(a)? )?dzx is positive, meaning it has no influence on min- 
imization and according to (9), the equation can be expressed as follows: 
1 
Imn(Z, F) = svec(Z)" (Pn(L)®n(T)" ® Tm) vec(Z) 
— H?(@,(T)? @ Ym) vec(Z) 


+ SF Unt1 8 (DY TmDY*))F, (32) 


where 
Tm = diag(ho,...,hm). 


The problem of optimal control in the discussion has now become a finite 
dimension optimization. We use the Lagrangian multipliers method to solve 
the ensuing optimization problem. Let us clarify 


J* (2, f) ~ J*(Z, FQ) = Jmm(Z, F) + AT (Mvec(Z) — NF-—Q), (33) 


where 
T 
A= [Ai Ag «+. A(mat)x(ne1) |, 


which shows the vector of Lagrange multipliers as A. The following equations 
lead to the following optimality conditions: 


OF* (2, f) _ 
Ovec(Z) 


OJ" (2, f) 
OF 


OJ" (2, f) 


7) ns 
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The Newton iterative approach or MATLAB software tools can be used to 
solve the aforementioned algebraic equation system. We can get the approx- 
imate solutions z(x,t) and f(x) from (18) and (19), respectively, by figuring 
out Z and F. 


6 Numerical examples 


This section gives two examples along with figures to illustrate how the rec- 
ommended technique may be implemented successfully and to show its po- 
tential. The results of the existing plan are analyzed and compared to the 
solution that was found analytically and method of [37]. The rand function 
is used by the MATLAB software to generate noisy data, and the value of h- 
for 0 <6 < 1 in the collocation points ety ae is calculated as follows: 


he =h+0.rand(size(h)), (34) 
m+1 
= aie _ pl2\3 
€ = [lhe — hile = (— 5 y |he — hI?)?. (35) 


For noisy data h-(x), the goal of the optimal control problem is to find 
functions f(a) and z(z,t) that minimize the following objective function and 
satisfy (1): 


a: 1 (oy ! 
min J(z,f) = 5 || 2(@,1) — he(@) IIi2p0,y +5 IF @) Mien 
i 2 of’. 2 
=5f le@1)-he()Pde +2 f lf (@)Par. (36) 
0 0 
We take the regularization parameter o = ¢’, and, in order to observe the 


convergence of the method described in numerical experiments, we calculate 
the approximate error resulting from the following equation: 


e(f) =||f-flln., (37) 


where f is the numerical approximation of the exact solution f in the collo- 


cation points {€}7t". 


Example 1. Consider the inverse problem with 0 = (0,1) x (0,1) [37] 
2i( Bt) — eax (2,1) = f(z), (x, t) € o, 
z(x,0) =0, x€ (0,1), (38) 


Z2(0,t) = 2,(1,t) = 0, te (0,1). 
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We attempt to approximate the heat source defined by 
f(z) = x? cos(r2). (39) 


Then with f given by (39), the forward problem presented by (38) has an 
analytical solution as follows: 


[oe) 


2(a,t) = y 


n=1 


Ls en (nn)*t 
5 Fn cos(n72), (40) 


(nm) 


where f,, is the Fourier coefficient as follows: 


fr 2 f f(x) cos(nrax)dx. (41) 
From (40), we have 
1 = en (nm? 
h(a) = 2(a,1) = s “Tame tn cos(nrax). (42) 


Table 1: Comparison of errors estimate obtained for functions f in Example 1 over a 
range of o values between the proposed method and [37] 


ei? e=10-* p= 10° 
Proposed method | 6.1387e —01 6.8907e —02 1.2413e — 02 
Method of [37] 2.51le — 01 4.8le — 02 3.26e — 02 


Table 1 analyzes the error behavior of the proposed method in here and 
the presented method in [37] by varying the value of o. 


Table 2: Errors estimate for the functions f and z in Example | over a range of o values 


ao=10-4 o=10° a = 10-8 ao=10-7 
Eroor(f) | 6.1387e—01 6.89076 —02 1.2413e—02 8.3175e — 03 
Eroor(z) | 5.5543e—02 6.0304e—03 6.0868e—04 6.3734e — 05 


The approximate solutions for the functions f and z are shown in Figure 
1. The approximation inaccuracy are shown in Table 2. Figure 2 depicts the 
convergence of the suggested approach. 
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Figure 1: Results of Example 1’s numerical solutions for functions f (left) and z (right) 


for a range of o values 
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Figure 2: Convergence of the numerical solutions of Example 1 for functions f (left) and 


z (right) for a range of collocation point values 


Example 2. Consider the inverse problem with 6 = (0,1) x (0,1) [28] 


2t(a,t) — Zea(x,t) = f(x), 


2(x,0) = sin(ra), 


z,(1,t) = —z,(0,t) = m(e~*™ * — 2), 


(x,t) Ee O, 
x€ (0,1), (43) 


te (0,1). 


We attempt to approximate the heat source defined by 
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f(x) = 2n* sin(rz). (44) 


Then with f given by (44), the forward problem presented by (43) has an 
analytical solution as follows: 


z(x,t) =—(e~” * — 2) sin(rz). (45) 
From (45), we have 


h(x) = z(a,1) = —(e~™ — 2) sin(ra). (46) 


Table 3: Values of errors for the functions f and z with different values of o in Example 
2 


g=10-° a =10-° og =10-7 a =10°8 
Eroor(f) | 4.3462e —00 1.2586e—00 7.4287e—01 6.8760e — 01 
Eroor(z) | 9.7406e —02 1.6259e—02 1.7443e—03 1.758le — 04 


Figure 3 shows the approximate solutions for the functions f and z. The 
approximation error is presented in Table 3. The convergence of the proposed 


method can be seen in Figure 4. 
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Figure 3: Behavior of the numerical solutions for the functions f (right) and z (left) at 
some different values of o in Example 2 
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Figure 4: Convergence of the numerical solutions for the functions f (right) and z (left) 
at some different values of collocation points in Example 2 


7 Conclusion 


When it comes to finding a regular and stable solution, inverse problems 
that are related to PDEs provide a significant computing challenge since it 
is very difficult to do so. The scope of this work is an investigation into an 
inverse space-dependent source issue for a heat equation. A shifted Legen- 
dre polynomial and an optimum control strategy were used in the process of 
creating a heat source. One of the most popular and efficient tools for resolv- 
ing computing problems is the Legendre polynomial. The shifted Legendre 
polynomials operational matrix was utilized to resolve this optimal control 
problem. By applying the suggested collocation method and using an oper- 
ational matrix, the issue was converted into a set of equations that can be 
solved using algebra. When utilizing this method to solve an inverse prob- 
lem, as demonstrated by the examples provided in the paper, a high level of 
precision was achieved in the solution. The method presented here was based 
on the optimal control problem and shifted Legendre polynomials. In future 
work, we will try to use machine learning techniques, including deep neural 
networks, to solve this problem. 
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